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ABSTRACT 

The feasibility of utilizing a gaseous core nuclear 
reactor to provide high enthalpy, high pressure gas flow 
for simulating atmospheric re-entry conditions was inves- 
tigated. The test facility uses a mixture of nitrogen and 
uranium in a closed cycle with no attempt to contain the 
uranium fuel within the core. The primary purpose of the 
facility is to provide high enthalpy, high shear flows 
for testing re-entry materials and shapes. 

Investigated in this study were the effects of the 
nitrogen-uranium mixture on reactor criticality, nuclear 
contamination of the test model, protection of the reactor 
core and nozzle structure from imposed heat loads and operating 
limitations of the test facility. 
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INTRODUCTION 



Many problems areas have been uncovered and investi- 
gated in our attempt to master space . One area is the 
re-entry of a ballistic, orbital or superorbital body 
into a planetery atmosphere and the delivery of an 
payload to the planet's surface. 

The kinetic and potential energy of the body must 
be dissipated to insure successful payload delivery. Atmos- 
pheric braking, with the transfer of energy from the high 

velocity re-entry body to the surrounding air, is often 

1 

the optimum entry system. 

The body must be protected from the high heat loads 
due to friction and the rapid compression heating of the air 
immediately surrounding the body. Ablative materials are 
used to both protect the body, by acting as an insulator, 
and to carry away the heat that is absorbed by the ejection 

. 2 

of mass of the ablative material. 

The ablative proq.ess, however, is very complex and 

attempts to analytically predict the behavior of new 

3 

materials and shapes have met with little success. Various 
attempts at simulating re-entry conditions have been made, 
but no existing facility can completely duplicate the 
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re-entry environment, although some aspects of the problem 

4,5 

can be simulated. 

Existing test facilities cannot duplicate the high 
enthalpy, high shear regime for a sufficiently long test 
time to be able to investigate this critical ablation region. 

This study will investigate the feasibility of a high 
enthalpy test facility which uses a fissionable gas reactor 
as a power source. The gas core reactor drives a closed 
cycle test system and provides the necessary energy for 
high enthalpy, high pressure gas flows of long duration. 

In establishing the feasibility of this system, the 
re-entry problem is first investigated to establish the 
power requirements of the system. 

Stagnation pressure and enthalpy are determined by 
the method outlined in Reference 1, from these the 
power required to simulate an atmospheric entry in a ground 
test facility is determined. 

Then the reactor itself is studied. The gaseous core 

reactor utilized is similar in operating principle to 

the gaseous core nuclear rocket engine proposed as the 

6 , 7,8 

next generation space booster. However, as the 

proposed system incorporates a closed operating cycle, 
containment of the nuclear fuel does not present an obstacle 
as it does in the gaseous core nuclear rocket. Critical 

size, fuel concentration and fuel-nitrogen ratio, where 

i 

nitrogen is the working fluid, are determined using the 



3 



diffusion approximation and one group neutron theory. 

The contamination problems’ associated with the utilization 
of a nuclear medium are also investigated. 

The combined heat loads to the reactor and nozzle 

10 

are examined. The heat loads in the nozzle are determined 
and homogenous transpiration cooling^'*' to protect the nozzle 
is examined. The use of a seeded carbon particle layer 

to protect the moderator from radiation heat loads is also 

• ^ 12 / 13 

examined . 

Finally the power requirements of auxiliary equipment 
to maintain steady state closed cycle operation are examined. 
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CHAPTER I 



RE-ENTRY CONDITIONS 



1,1 Pressure and Enthalpy 

To test the materials used in a re-entry vehicle it 
is necessary to know the velocity, the pressure and the 
heating involved in the re-entry. 

As a typical example consider a vehicle returning 
to earth with a superorbital velocity. An estimate of 
the velocity can be obtained by assuming it equal to the 
velocity needed to escape from a nonrotating earth, with 
no drag forces. The energy per unit mass to escape from 
the earth, acting as a central force field is 




Ik. 



(l.D 



and the velocity to escape is 



V<£ 




( 1 . 2 ) 



when k is the earth gravitational constant, 4 mo** '"VseS 
and R> is the earth's radius , G.4. > io & *r\ . , 



the escape 
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velocity from the earth's surface is 1.12 x 1C) m/sec. 



This is the minimum for earth escape and is the mini- 
mum at which a superorbital body will approach the sensible 
atmosphere, commencing re-entry. 

Expressions for stagnation pressure and enthalpy, at 
the point of maximum pressure and also at- the point of 
maximum heating rate, may be found by using the velocity 
altitude history given in Reference l for an exponential 
atmosphere, 

po « e ^ 



v = v c e 



i 



- Cp A 

2 Sm ©c 



ys 



} 



(1-3) 



where V«. is the initial entry velocity, Co is the coeffi- 
cient of drag of the body, the mass of the body, ©e 

the entry angle and $( and y3 are constants defined by 



f • f. 



e 



(1.4) 



The stagnation pressure on the body is approximated 
by Newtonian theory with the body angle considered small 
and the ambient pressure neglected. The stagnation pres- 
sure at any point in the entry is 



f© * ^> 0 Q ^ 



-CoK 






Po < \ 

) 



(1-5) 



cA. Po 



* o 



For maximum stagnation pressure, set 
obtain 



and 
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f- -C l e^('-/3)e CiejS \ C.e^e Cie ^ C,eMft) *0 

where 

C, = f> b ( ^X/e 1 

^2 - _ C P A «_£o 

m-Au\Oe. 



( 1 . 6 ) 



C v and Cz. are fixed by the atmospheric model 
ballistic parameter, an< ^ initial entry velocity; 
and remain constant during re-entry. 



The altitude for maximum stagnation pressure during 
re-entry is given by 

^ d. 7 ) 

Since the velocity for maximum stagnation pressure is 
given by 

V = V« e * (1-8) 



and associated density is 



p _ yy>S)A 6e A 

‘ C D A ' 



the maximum stagnation pressure is 

P _ m S\r» 0 « 



( 1 . 9 ) 






Co k 



( 1 . 10 ) 
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and the stagnation enthalpy is 



Mo - la ♦ 



(l.lla) 



or since the static enthaply is very much less than dynamic 
enthalpy, 

i ,2. 

(1.11b) 



u S 

Ho 2e 



Similarly, expressions are obtained for stagnation pressure 
and enthalpy at the point of maximum rate of convection 
heating of the stagnation point, where stagnation point 
heating rate is given by^" 

S'- C */F v 5 



( 1 . 12 ) 



where is the stagnation point coefficient of heat trans- 
fer, a function of body configuration, <T* is the radius of 

curvature at the stagnation point, and ^ and V are functions 

of altitude. To obtain the point of maximum stagnation point 
heating rate, the previously derived expressions for 
and V are substituted in the above. 

<J t ’ 

For the maximum set — ( dLt - y = 0, obtaining 



SO 

which occurs at an altitude of 






' ©e 


A 




(i.i-3) 


do k 


3e<r 




tt 




(1.14) 


W Q (2 


J 
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and velocity of 



_ 7c, 



v ' v« e 



(1.15) 



The associated stagnation pressure is 




r~r> S\ A ©e 

C& & 




(1.16) 



and stagnation enthalpy is given by 




(1.17) 



Figure i shows the extremes of stagnation pressure 
verses stagnation enthalpy which are encountered during 
an atmospheric re-entry. The ballistic parameter, 6 c 

is a free parameter in determining the envelope. Two 

/ 4 

envelopes are shown, one for ve- 1.12* 10 ta/sec. t corresponding 
to velocity entering the sensible atmosphere equal to 
escape velocity. The other is for Ve = O.© * lO* vn /sec., 
corresponding roughly to lOOrt.rt. orbital velocity. It 
can be readily seen in Figure 1 that: 

i) stagnation enthalpy at maximum heating rate and 
maximum stagnation pressure is independant of 
the ballistic parameter, but increases with in- 
creasing initial entry velocity. 

ii) stagnation pressure increases directly as entry 
velocity and directly with ballistic parameter. 

1.2 Power Requirements for Simulation 



The stagnation pressure on the test model is fixed 
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by re-entry conditions but the chamber pressure is limited 
by structural considerations. The maximum chamber pressure 
is set at 1000 atmospheres for the analysis. It is desir- 
able to test at the lowest possible Mach numbers consistent 
with representing hypersonic pressure and heating rate dis- 
tributions. For this reason, and to minimize nozzle losses 
and cooling requirements, a test Mach number, K-t , of 3.0 
is used in this analysis. 



For isentropic frozen flow through the nozzle, the 
power required per unit throat area is 




PoMo 

\/2RT* 




(1.18) 



If Equations 1.16 and 1.17 are used for "fo and 
^\ 0 respectively, the power required per unit throat area 
to simulate the conditions at the point of maximum heating 
rate is 



£_ , J V«/*- 

A* &e 2/3 / 



(1.19) 



These power requirements are plotted in Figure 2 against 
initial entry velocity with the ballistic parameter as a 
free parameter. Reference 14 was used to obtain the charac- 
teristics of nitrogen, which was taken as the test gas. 

Limiting the chamber pressure to 1000 atmospheres 
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established one limit on the re-entry profiles that can 
be simulated. Another is- set by limiting the chamber 
temperature to 10, 000° K, this upper limit is set after 
a consideration of the heat transfer loads experienced 
in the various test facility components. 

1.3 Test Time 



The main purpose of this facility is to examine the 

behavior of re-entry materials and shapes under the high 

heat and shear loads of an atmospheric re-entry. Such 

materials are ablative in nature, that is, they absorb 

heat by increasing in temperature and changing in chemical 

or physical state. The heat is carried away from the 

surface by a loss of mass. The departing mass also blocks 

part of the convective heat transfer to the remaining mate- 
2 

rial. The tunnel running times are set by the time required 

for the rate of ablation to reach equilibrium. This time 

15 

is a function of the ablative material and for typical 

2 

ablative materials , such as, graphite and glass reinforced 
phenolics, the ablative equilibium times is extended 
to 200 seconds 

1.4 Existing Facilities 

Shock tubes are used to generate higher temperatures 

and pressures to simulate re-entry conditions. The Cornell 

16 

Aero Laboratory six foot hypersonic tunnel can generate 
stagnation temperatures and pressures up to 4200°K and 
2000 atmospheres respectively. The test time for these 
pressures however is only 4 micro-seconds. This tunnel can 
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be upgraded to 10,000°K and 2000 atmospheres but still has 
a run time far too short to produce the time history of 
high temperatures and pressures at which the ablative 
materials fail. 

The test time for a shock tube may be increased by 
using several tubes fired in a Gatling gun fashion. This 
however only increases run time to seconds, developing 
a jet power of 25,000 HP , about four orders of magnitude 
too small for re-entry requirements. 

High enthalpy and pressure flow is also possible 
utilizing a hot shot tunnel powered by an electrical 
discharge arc. This type of tunnel requires large amounts 
of stored electrical energy and is presently limited to 
a few hundredths of a second operation.^ 

Although these facilities can approach the desired 
levels of enthalpy and pressure, they cannot maintain these 
conditions for a sufficient lenght of time. They are severely 
power limited. As an example, the power required, from 
Figure 2, to simulate the re-entry of a vehicle whose ball- 
istic parameter is 2000 grams/ cm and initial entry 

2 

velocity is 8000 meters/second , is 30 megawatts per cm of 
throat area. This is about ten percent of the total output 
of the Grand Coulee Dam for each square centimeter of nozzle 
throat area. 

Plasma arc tunnels provide high enthalpy but low 
pressure flows and are currently being used to test re-entry 
materials. These facilities are also power limited and operate 
at low stagnation pressures and use argon as a working fluid 
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to achieve long duration tests with comparatively low 
power inputs.^ 

Because of these high power requirements, a gaseous 

core reactor is suggested to drive the test facility, the 

gaseous core because its operating temperature is limited 

only by heat loads to the containing structure and can 

o 

operate at the desired temperature, 10,000 K. 
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CHAPTER II 
TEST FACILITY 

2.1 Gaseous Core Reactor Requirement 

The power required to produce long duration high 
pressure ,high enthalpy gas flows appears to be technically 
achievable only with nuclear power. Conventional combus- 
tion is limited by the recombination of combustion products, 
which absorb energy and result in an equlibrium tempera- 
ture far below the required stagnation temperature. In 
atomic reactors, solid and liquid core operating tempera- 
tures have upper limits, the former is the melting point 
of the nuclear core and, or, structure and the later is 
boiling point of the liquidized fuel, and, or , melt ing point 
of the structure. 

A gaseous core nuclear reactor is capable of steady 
state operation at very high temperatures and pressures. 
The operating temperature of a nuclear reactor utilizing 
a fissionable gas core is limited only by the heat load 
that the hot gas imposes on the container. These heat loads 
can be controlled to some degree by a combination of con- 
vective and transpiration cooling and by injection of a 
thermal radiation absorber to block radiation heat transfer 
to the wall. The gaseous core reactor in this study is 
capable of providing the energy necessary to run a test 
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at the desired stagnation pressures and enthalpy, that is, 
7,100 megawatts. 

2.2 Gas Cycle 

To simulate the effects of atmospheric heating, nitrogen 

is used as the working medium. Uranium 235, initially a 

fine power, is the nuclear fuel. Biological hazards 

and economic necessity make it mandatory that the uranium 

fuel be contained within the core of a gas core nuclear 
18 

rocket. As the proposed system is a closed system, these 
considerations do not apply. This is not to say, a priori, 
that fuel containment in the core proper is not required. 
Presence of uranium in the nitrogen flow complicates the 
analysis and interpretation of experimental data. If the 
ratio of mass flow of uranium to nitrogen is limited to 
a maximum value of ten percent, or about one percent mole 
ratio, it should permit the assumption that the flow in 
the nozzle and test section behave as a pure nitrogen flow. 
This limit is imposed as the. maximum allowable uranium 
content for thermodynamically simulating the effects of 
atmospheric entry. 

The presence of fissionable gas in the closed cycle 

causes induced radioactivity of the system components. 

State of the art shielding and radioactive materials handling 

q 

techniques are assumed sufficient to provide necessary 
biological protection. As this facility would be very 
expensive to construct, it is necessary that the induced 
radioactivity of the model be within biological limits 
for examination of the model within a short time after 



the test and associated radioactive exposure. 



The radioactive damage or effects to the model must 
be slight enough so that these associated effects do not 
mask or override the aerodynamic effects which are being 
investigated . 

These nuclear effects on the model are investigated 
to determine their influence on the design of the over-all 
system, and to determine if containment, greater than that 
implied by the thermodynamic limit of ten percent .uranium, 
is required.’ 

Chamber pressure and temperature are set at 1000 
atmospheres and 10,000°K. A portion of the heat load to 
the moderator is absorbed by the flow of nitrogen entering 
the core, the remainder is removed by an external water flow 
enroute to another system component, the water diffusor. A 
seeded layer of carbon is used to provide a radiation barrier 
which protects the moderator from radiation heat loads. 

The thermal energy of the hot gas in the reactor core 
is converted to kinetic energy in a convergent -divergent 
nozzle. The nozzle is cooled by conventional backside 
cooling and transpiration cooling, both utilizing a secondary 
flow of high pressure, but low enthalpy nitrogen. 

A schematic of the test facility is given in Figure 3. 

A temperature - enthalpy diagram of the nitrogen cycle is 
given in Figure 4. 
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After the test section, the secondary nitrogen flow 

is vented into the main flow, reducing the temperature 

slightly. The stagnation temperatures are still in the 
. . . o 

vicinity of 10,000 K and prevent the use of conventional 
turbo-expansion devices or heat exchangers at this point. 

A water spray diffusor is utilized to slow the flow 

and lower the stagnation pressure and temperature of flow 
o 

to 1000 K and 100 atmospheres. The water spray is sufficient 
to cool the walls of the diffusor. 

A centrifugal separator is used to separate the 
condensed uranium particles from the flow of superheated 
steam and nitrogen leaving the water diffusor. 

At this point the flow is basically free of uranium 
and at a temperature and pressure at which conventional 
turbo-machinery can be employed. The flow has 2,500 megawatts 
of available power, if expanded only to the saturation point 
of the water. After this energy is removed from the stream, 
the water is condensed and separated. The nitrogen, at 
0.24 atmospheres and 370°K,is cooled further in a conven- 
tional heat exchanger and pumped back to 1000 atmospheres 
using two multistage centrifugal flow compressors with 
interpump cooling to maintain temperatures below 1000°K. 

The power required for the pump work, 880 megawatts, is 
provided by the aforementioned turbine . 

The pressure of the water is raised to 120 atmospheres 
in another multistage centrifugal pump, requiring 10 mega- 
watts of the above turbine work. 
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The temperature - enthalpy diagram of the water cycle 
is given in Figure 5. 

2.3 Component Interface Problems 

To be successful as a system, the components of the 
facility must be mutually compatible. In the following 
sections, the operating characteristics of the reactor 
and the nozzle, which are the main components of the 
system, are examined. The reactor is examined to determine 
the effects of various operating pressures and fuel-nitrogen 
concentration on the critical size and fuel concentration. 

Contamination of the model is examined to determine 
if it imposes a constraint on the fuel-nitrogen mixture 
in the core. 

Heat transfer in the reactor and nozzle are also 
examined to determine limits on operating temperature, 
power and component size. 
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CHAPTER III 

GASEOUS CORE NEUTRONICS 

3 . 1 General 

Reference 6 shows that the gaseous core nuclear 
reactor is theoretically feasible. It shows that one 
can obtain a reasonable estimate of reactor critical 
radius and critical fuel concentration by considering 
the reactor core to be a homogeneous mixture of fuel 
and moderator or a reflected sphere of fuel alone. The 
effects of containing the fuel with a gaseous annulus 
of hydrogen are also investigated as were the effects of 
multiple cavities. Reference 6 quantitatively shows that 
a gaseous core reactor can be designed with either simple 
or multiple cavity geometry. 

In this study, the reactor will be considered to be 
a homogeneous mixture of pure uranium 235 and Nj. with 
a bare spherical geometry. The effects of reflecting 
the above bare sphere are also considered. One group 
neutron theory is utilized, that is, the neutrons are 
considered to be in thermal equilibrium with the reflector- 
moderator . 

No attempt is made to separate the fuel from the 
working fluid, nitrogen. Criticality calculations , 
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are made with the ratio of number density of nitrogen 
atoms to number density o’f fuel atoms, / I j 



as a parameter. This gives one measure of contamination 
of the flow in the model test section. From the fluid 
dynamic point of view, the flow is treated as behaving 
as pure nitrogen, and the number density ratio gives a 
measure of how valid this approximation is. The model 
is also subjected to radioactive contamination, the 
number density ratio gives a partial measure of this effect; 
partial in that the uranium itself in the flow is not 
the only source of model contamination. 

3.2 Criticality Condition . 

As a first approximation, the reactor can be considered 
to be a bare homogenous mixture of pure uranium 235 and 
nitrogen. Criticality is determined by considering all 
of the neutrons in the core to be thermal, that is, there 
is no loss of neutrons in the process of slowing down 
from the high energy state of production in the fission 
process, to capture at a lower thermal energy state. For 
steady state operation, the conservation of neutrons requires 
that the rate of production of neutrons in a control volume 
be equal to the rate of loss of neutrons due to absorption 
and leakage from the volume. 

The rate of production of neutrons per unit volume 
is given by 




Production 



3.1 
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where ^ is the average number of neutrons 
produced per thermal fission 
0 is the neutron flux 

and is the macroscopic fission cross section 

of the core material. 

The neutron leakage is given by the diffusion approxi- 
mation as 

Leakage = - D * 3 *$ (3.2) 

where V 2 is the La Placian operator 

and D is the diffusion coefficient, 

1/3 of t-he transport mean free path 
in the reactor 

The absorption term is similar to the fission term. 

Absorption (3.3) 

here is the macroscopic absorption cross section 

in the core. 

Equations 3 . 1, 3 .2 and 3.3 are combined and 
the conservation of neutrons is expressed in differential 
form by 

Dvy V -£o.) (/> = O (3.4) 

or if the material buckling, , is defined as 

8~ - 0(2*-^)/ o 



(3.5) 
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the conservation of neutrons becomes 



+ B* $ = O (3.6) 

For a spherical reactor, with the boundary conditions 
that the flux goes to zero at the radius of the core and 
the gradient of the flux is zero at the center of the 
sphere. Equation 3.6 can be solved for the geometric 
buckling of the reactor. 



a 1 = 



(3.7) 



The criticality condition is the matching of the 
geometric buckling which satisfies Equation 3.6 with 
the material buckling of Equation 3.5. The diffusion 

length is defined as 



L 1 = 

and the infinite multiplication factor as 

= V 3 



(3.8) 



With this notation, the criticality requirement is 

1 _ fcc -1 

L* 



WrY 



(3.9) 



(3.10) 



The criticality requirement for the system is that 
the geometric buckling be equal to the material buckling. 
This means that the rate of neutron creation due to fissions 
is equal to the rate of destruction of neutrons due to 
absorption, that is, fission and nonfission capture, and 
leakage from the critical region. 



For a nitrogen and pure uranium 235 system, the number 

18 

of neutrons released per thermal fission of uranium 235 is 
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v>= 2.4fc jp**-?,* 



»-»s 



The interaction cross sections per atom for uranium, 
sub-subscript u. , and nitrogen, sub-subscript , are 
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G*«. = 10 barns 
S'** = 10 barns 

where the subscripts 4* denotes fission, 

gl. denotes absorption, 

& denotes scatter. 

With these cross sections, the criticality equation 
can be used to determine the number density ratio of 
nitrogen atoms to uranium atoms required for criticality 
as a function of the core radius and the number density 
of uranium 235 atoms. 



= 580 barns 
<r = 680 barns 
=1.88 barns 



Mu 



- 199 rywl 



47 * 

t-SAS X ID 

R* Mo 1 



(3.11) 



Using this relation the critical radius is found as 
a function of uranium 235 concentration with >4** /t4o 
as a parameter. The results are given in Figure 6 with 
the mass flow ratio of nitrogen to uranium shown in lieu 
of the number density ratio. 



Examination of the criticality equation above shows 
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that as Mn / Mo decreases below 199, ft hi must increase. 
Therefore, Km /nJj = 199 gives the minimum reactor size 
for a given concentration, and the associated mass flow 
ratio is 15.4. This is due to the scattering contribu- 
tions of the nitrogen atoms in the core reducing the 
transport mean free path of the neutrons . This is most 
readily seen if the criticality condition is expressed by 

O*- _ TT 2 3(1- 4046^ 

K KlJ 1 l V> V, - 03L* - if- K <T*o ^ ’VU- ) (3.12) 

T 

where the macroscopic cross sections are replaced by 
the number density times the microscopic cross section. 

The denominator has a maximum value when Km / is 

199. For higher or lower number density ratios, the 
critical radius increases. 



As used here ft. is not the actual radius of the 
homogenous sphere, but rather that radius plus the extrapolation 
distance. The extrapolation distance is the distance 
at which the neutron flux, which has a finite value at 
the physical boundary of the sphere, goes to zero. Reference 
9 gives the extrapolation distance as 

J=2.rsD (3.i3) 

The physical radius of the bare sphere is 

r » R- z.13 0 (3 - 14) 

This critical size and uranium concentration so 
obtained is only an approximation and is valid only- in 
determining the effect of nitrogen poisoning of the reactor 
on critical size and critical uranium concentration. The 
actual reactor will be operating at elevated temperatures 
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which will effect the interaction cross section, which 
generally vary inversely with neutron energy. A higher 
equilibruim temperature means higher neutron energies 
and, in general, lower interaction cross sections. Also 
the actual reactor will not be bare but must be enclosed 
in some way. 

The reactor core is in fact contained in a pressure 
vessel, for structural reasons the walls of which must be 
maintained at a low temperature relative to the reactor 
core. Lining the pressure vessel with beryllium, a 
moderator, which must also be cooled for structural reasons, 
provides a means of decreasing leakage from the core and 
a means of providing low energy neutrons to the core . 

Nitrogen, being heavier than beryllium, does not 

i 

absorb as much of the neutrons energy in a collision, 
therefore more nitrogen collisions are necessary to thermalize 
a neutron. Also nitrogen, even at very high pressures, is 
still* very much less dense than the solid beryllium moderator 
and at the higher core temperatures, the scattering cross 
section for nitrogen is less than that for the cold beryllium 
moderator. It is reasonable to assume then, that the 
neutrons are in thermal equilibruim with the cool reflector 
and not with the hotter nitrogen-uranium mixture. 

The effects of the reflector-moderator on core radius 
can be found from a one velocity, two region consideration. 

As before in the core, the conservation of neutrons is 
given by Equation 3.6. 
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In the reflector, there is no fissionable material 
and the conservation of neutrons gives 

Dr ~ = O ( 3 . 15 ) 

where subscripts r and c stand for reflector and core 
respectively - 



With the inverse of the square of the reflector 
diffusion length, L. r , defined as 




(3.16) 



the neutron conservation equations in the core and 
reflector is obtained by dividing Equation 3.15 by Dr / 



< 7 * (fic. + Be cf><_ =0 (3.17) 

and 

V* - k. r = O (3.18) 



This set of differential equations can be solved 
for a criticality condition with the boundary conditions 

O < ^ c (r = o') < co 

•V $cC r ~R.*') c 

= D - V ^L*. (3 - 19c) 



(3.19a) 

(3.19b) 



and 



^r(Rr + dr + S') 



O 



(3 . 19d) 
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where d>- is the reflector thickness 

and R,. is the critical radius of the reflected core. 

The solution of Equation 3.17 and 3.18 with 
boundary conditions of 3.19 gives an expression for 
the criticality of a reflected spherical reactor in the 
transcendental form, 

Q.(k.R v- CoWi ( c\«- +3 ) t \ ) - Dc ( l "Rr Be Co\ Be R<0 (3.20) 



The reflected core radius can be expressed in terms 
of the bare core radius, r , and the reflector saving, S, 
as 



R,- * r-S 



(3.21) 



This value is substituted into Equation 3.20 and 

a transcendental expression for reflector saving is 
obtained. 



col £S 



Per V (dr - 0 - ) 



(3.22) 



Equation 3.22 is solved graphically for the reflector 
saving at each known r, then the reflected core radius, R r< is 
obtained by Equation 3.21 



The resulting critical radii are plotted in Figure .7 
against critical uranium concentrations for a 100cm 
beryllium reflector. Shown is the reflected core radius, 
with the mass flow ratio of nitrogen to uranium as a free 
parameter. Also shown are lines of constant pressure for 
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an operating temperature of 10,000 K. 

Heat transfer in the reactor core, coupled with 
practical limits on the capacities of the system's pumps 
and turbines, lead to the choice of a minimum core radius. 

Structural consideration dictates a low operating 
pressure, but this is not compatible with the desire 
to simulate the high stagnation pressures of a re-entry. 
Also, the problems associated with heat transfer in the 
core require high mass flows and low temperatures. For 
these reasons the operating pressure is set at the maximum 
possible with projected structural technology, that is, 

1000 atmospheres. 

While it is also desired that the mass flow rate of 
nitrogen to uranium be maximized so that the test facility 
simulate atmospheric conditions, the requirements for a 
minimum core size dominate in this case. It is also noted 
that the critical core radius increases rapidly with a 
slight increase in mass flow ratio from the minimum of 15.4 

For these reasons the core radius is set at the minimum 
possible with a pressure limit of 1000 atmospheres. This 
radius is 170 centimeters, occuring at the mass flow ratio 
of 15.4 
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CHAPTER IV 

MODEL NUCLEAR CONTAMINATION 



4.1 Radiation Damage 

In the gaseous core nuclear reactor used as a power 
source for this system, there is no provision made in 
the core itself for the separation of fuel and fission 
fragments from the working fluid. This makes the design 
of the reactor much more simple then if some hydrodynamic 
means of separation of fuel from the working fluid is 
utilized. As the system is closed there is no problem 
of venting radioactive material to the atmosphere , relieving 
one of the major constraints that faces the use of a gas 
core nuclear rocket. 

The presence of radioactive material in the fluid 
stream causes contamination of the model being tested. 

This contamination is examined to determine if the structural 
damage to the model will signicantly influence or mask 
the aerodynamic and thermal effects being tested. As the 
model must be examined after the test, the extent of radio- 
active contamination as a biological hazard also is examined. 

Radiation damage is caused by neutron bombardment 

20 

and interaction, alpha and beta bombardment and gamma 
radiation . 
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The radiations may be divided into two groups. The 
light group consisting of betas, gammas and other electro- 
magnetic radiations, and the heavy consisting of neutrons, 
fission fragments, accelerated ions, and alphas. As a 

general rule the damage from the heavy particles is more 

2 0 

severe that that from the light particles. The heavy 
particles displace the atoms of a solid crystal lattice 
from their normal position. In general this will result 
in an increased ultimate strength, a decrease in elongation 
and a reduction in area of the material, also the micro- 
hardness of the material increases and there is a decrease 
in impact strength. These effects are minimized if the 
neutron bombardment occurs at elevated temperature, as 
is the case in this test facility. ■ 

Graphite, a prime candidate for a re-entry material, 
incurs crystallite damage during irradiation. Hardness 
and strength are increased while thermal and electrical 
conductivities decrease. Graphite shapes change in gross 
dimension during radiation. All of these effects decrease 
in magnitude when the temperature of the graphite during 
radiation is increased, with almost no effect when graphite 
temperatures are above 1000°C. Also the integrated flux 

levels at which these effects begin to be noticable is 

>9 . 2 ( 

high, of the order of 10 neutrons per square centimeter. 

As this test facility will provide very high model tempera- 

15 

tures and neutron fluxes at the model of the order of 10 
neutrons/ cm 1 sec, it is expected that irradiation damage 
will have negligible effect on the test results. 
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4.2 Biological Hazards 



In determining the biological radiation hazard due 
to model contamination the following assumptions are 
made : 

i. radioactive particles, such as, fuel, fission 
fragments do not adhere to the model upon 
impact , 

ii. contamination is due to neutron-model 
interaction causing a transmutation to 
a possibly unstable isotope, 

iii . contamination is acceptable if the model 

can be handled with minimum protection within 
a short time of the test run, 

iv. acceptable radiation levels from the model 

_3 

are 7.5 x 10 rem/hour as recommended by the 

21 

US National Committee on Radiation Protection. 

4 ,2a contamination Environment 

Neutron density at the model has three sources; 
direct streaming from the reactor core, neutron flux 
from the fissionable mass in the test section, and secondary 
emission from fission fragments. 

Direct streaming of neutrons from the core can be 
measured by noting the attenuation of the neutron flux 
by the gases in the nozzle between the model and the core. 
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The differential change in neutron flux due to absorption 
by the gas is 

af = 02- t 4 - 1 ) 

The macroscopic absorption cross section of the gas 
varies with position in the nozzle as the density of 
the gas decreases. Assuming this variation is linear 
with axial nozzle position, % , and that the density 
goes to zero at the model. Equation 4.1 becomes 

sg- = $ ( i-t:) d% (4 - 2 > 

When integrated and evaluated at the model, the 
flux at the model is given by 

, ~ ^/2- A a. o 

^ e ( 4 . 3 ) 

where L is the axial position of the model and Xo. 
is the absorption mean free path 

The axial position of the model, L» ,is limited 
by heat transfer considerations in the nozzle and inlet 
and must be minimized. Then Li is the order of the 
required nozzle length to expand the flow to Mach 3, 
approximately 30 cm, for the throat size set in Section 6.1. 

For the chamber condition of this study the mean 
free path for absorption in the chamber is about five 
times this value, then the free streaming flux from the 
core to the model is approximated by 



& = 0.9 & 



(4.4) 
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The neutron flux from the mass of fissionable material 
in the test section is found by considering the change in 
flux as the flow goes from the critical region of the 
core to the subcritical test section. 



The time of flight from the core to the test section 
is of the order of the length divided by the speed of 
sound at the nozzle throat, 

t f = L*/C* (4-5> 



If the test section is considered to behave as 

an infinite cylinder with constant cross section, the 

critical buckling is a function of the geometry and is 

9 

related to the test section radius by 




2.405 \ 2 

K* j 



(4.6) 



is 



The effective multiplication factor in the test region 




(4.7) 



where VCoo is defined by Equation 3.9 



With the assumption that the number of neutron generations 
is given by the time of flight, tf t divided by the average 
neutron lifetime, Acl/iXa / the neutron flux at the model, due 
to the reactivity of the fissionable gas in the test section, 
is given by 





e 



(4.8) 
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Equations 4.5 through 4.3 are evaluated for 

the conditions in this facility, 

_ A 

tf - 5* /O seconds 
keg m - 00Z 
5oo 

to obtain the flux at the model due to the reactivity of 



jg& zL. , 



the test section / 



-»5o 



“ & € (4.9) 

As is very much smaller than 9^5 it 

is considered negligible in this analysis. 



Secondary emission of neutrons by fission fragments 
in the vicinity of the model will also contribute to the 
neutrons bombardment of the model. Given a fission fragment, 
A , with a half life, ty z A , and a yield fraction 
per fission of # from the definition of half life, the 

radioactive decay constant t t is 



^ m J/rt.2. 

*\ "w 

The rate of creation of atoms of fragment A is 



(4.10) 



(4.11) 

where £y^is the number of fissions per unit volume per 
unit time. 

The rate of decay is 



Thus in the 
the rate 



4* ; -Xa A, 

dl 



reactor, fission products. 



A 



, build up at 



(4.12) 




A 



(4.13) 
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and concentration of species A at time i is 

A M - t (m» 2, <k - {r[» irf. e' >Ai ) (4, 14) 

With the initial condition that the initial concentration 

of species A is zero, the concentration of fragments leaving 

the core is given by A (to) , where t© is the dwell time 

of fluid in the core, 

L , (4.15) 

At t© =1 second, the concentration of A leaving the 
reactor. is given by 

( 1 - 0 ^ ) (4.16) 

Aa 

The rate of decay in the nozzle is given by Equation 4.12 
evaluating this at the test section with a time of flight, 
t$ , from the core, the rate of decay of fragment A 
at the model is 



* - Xa Ik tn e ' A ") e (4.17) 

If it is assumed that one neutron is emitted per 
decay, the sum over all fission fragments of the decay 
rates given by Equation 4.17 is the rate of creation 
of neutrons at the model, 

2 ft* O' 6 "* 1 ) (4.i8) 

A 

The absorption mean free path in the test section 
is about fifty times larger than the test section radius, 
therefore all neutrons freed in the test section will hit 
the model or walls. If all neutrons hitting the walls 
are absorbed, the ratio of neutrons hitting model to those 
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hitting the wall is given by the ratio of model radius to 



model blockage of the test section assume this ratio is 0.1 

Then the delayed neutrons strike a unit area of 
the model at a rate 



The fission products yield are given in Reference 22. 

An examination of the distribution of fission products 

shows that for a mass number less than 80 or greater than 

160, U* is less than one part in ten thousand and the 

contribution to $*a is negligible regardless of the half 

life. Also, if is greater than one minute, is less 

than .01 seconds \ this combined with a maximum YJ* of 

ten percent from Reference 22, shows that the contribution 

of any element with a half life much greater than one minute 

may be ignored in approximating the neutron flux due to 

secondary emission from fission fragments. References 20 and 

22 were examined for fission products with yield fractions 
-4 

greater than 10 and half lives less than one minute. The 
isotopes satisfying these criteria are given in Table 1. With 

these values, the rate of creation of neutrons at the model is, 

-2 

by Equation 4.18, 1.72 x 10 neutrons per second and the 
neutron flux at the model due to secondary emission of neutrons 
is 



test section radius, Rm / <8.* 



For negligible 




(4.19) 



= 2.3 x io" 4 <$ 0 



(4.20) 



It is noted that this method gives a conservatively 
high estimate of flux at the model due to assuming that 
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the total fission yield for a given atomic number indicated 
in Reference 23 is the element with the atomic number with a 
short half-life, and further that these radionuclides all 
emit a thermal neutron upon decay. 

There are two unconservative assumptions in this 
development, that is that the initial condition is the 
absence of fission fragments at the beginning of the reactor 
dwell time, and that only short lived neutron emitters are 
involved. In operating as a closed cycle , the supply of 
nitrogen and uranium will be contaminated by the fission 
fragments generated during previous residence in the reactor 
core, increasing the initial concentration of fragments and 
allowing build up of fission fragments with a long half-life. 

Due to the shear magnitude of the power involved in 
this facility, it is expected that the time for a given mass 
of fluid to complete a cycle will be long enough to allow 
decay in the concentration of the prompt neutron emitters 
and that provision can be made for monitoring and removal 
of long lifed fission fragments as a part of the fuel handling 
process . 

With these assumptions, the flux due to secondary 
emission from fission fragments is several orders of magnitude 
smaller than the direct streaming flux, and the former is 
neglected. 

The flux striking the model is then approximated by 
the direct streaming neutron flux of Equation 4.4 . 
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The average flux in the gaseous core is set by the 
reactor size, the power required, the energy per fission, 6, 
and the macroscopic cross section for fission of the fuel in 
the reactor core, 

P/Vo. 

< 4 - 211 

At a power level of 7100 megawatts with a 1.7 m core, 

15 2 

the average flux in the gas core is 3.26 x 10 neutron/cm sec. 

Equation 4.4 gives the flux at the model -for these 

. . 15 2 

conditions, 2.93 x 10 neutrons/ cm sec. 

4.2b Model Contamination 




The radiation hazard from the model is assumed to consist 
of secondary emission from neutron flux induced radionuclides 
on the model. The governing equation is of the same form 
as the build up of radionuclides in the core, then decaying 
after exposure. The build up of species A on the model during 
a test run of time, , is 



Mi,) 




(4.22) 



where is the neutron macroscopic absorption cross section 

in the reaction ft * 0 n‘ — A and B is one of the original 
materials of the model or a transmuted species. The rate of 
decay of A after the test is 



dA 

dt 



cMi-e^) e' w 



(4.23) 



where t is the time since exposure to the neutron. After 
the decay rate for a species has been determined, the type 
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of decay must be examined to determine the mode of decay, that 
is alpha, beta, <^amma or neutron emission and the energy of the 
emitted particle. The particle type and energy set the maximum 
decay rate that a man can be exposed to and not exceed the 
limits of .006 rem/day. It is noted that if more than one 
species is formed or more than one mode of decay is present, 
the sum of the exposures will determine the allowable decay 
rates of each species or mode. 

As graphite is the prime candidate for material to be 



tested, it will be examined. Carbon undergoes the reactions 







+ 


*n‘ = 




stable • 


(4.24a) 




4C' J 


+ 


on' = 


.C U tvx 


=5570 years 


(4.24b) 




c' 4 

6 C 


+ 


®n' = 


/ 3 ~ + T N* 


stable 


(4.24c) 


for 






8.55 x 


22 3 

10 atom/cm 







r- _ -27 2 

= 3.3 x 10 cm 

-27 2 

<r* = .9 x 10 cm 



t T = 200 seconds 

(zt t ) (4.25) 

The initial decay rate for the model is, by Equation 4.23 

Q ^ xiO T 

dll: Set - 

With no absorption in the model, the surface beta flux 

for a cylinderical model is conservatively given by 

d Cia % clCu R* 

dlt Am Ott 2 . 

For a carbon model with a radius of 1.25 centimeters 

d . 5 27 * \ 0 " T L^l*s*l/L 

dt sec 



(4.26) 
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Reference 20 gives the exposure of betas of various 
energies which produces an allowable dose of .006 roentegen 
equivalent man per eight hour day. The exposure rate 
above is on the safe side, therefore a graphite model may 
be handled immediately after a run at 7,100 Megawatts 

for 200 seconds. 

Correlations between nuclear decay products and allowable 

9,20,23 

dose rates are readily available in the literature 

A very slight protective layer, such as paper or cloth, will 

provide protection against alpha particles, which do not 

9,20 

present an external radiation hazard. 

It is also probable that high percentages of hydrogen 2 
will be present in an ablative material, and the decay of 
hydrogen isotopes yields a low energy beta particle in 
concentrations, determined by Equation 4.22 ; five orders of 
magnitude below carbon and therefore, as a graphite model, 
also may be handled immediately after a test. 

This analysis can not possibly determine if each and 
every material will be safe to handle after a test. The 
user will have to examine the reaction of each model 
individually to determine the contamination from each material 
and the time history of resulting total model contamination, 
as given by Equation 4.22 summed over all model components. 

The resulting model activity can be compared to allowable 
dose rates to determine if and when the model can be handled. 

It is noted by examination of Equations 4.4, and 4.22 
that the decay rate at the model is directly proportional to 
the flux rate in the reactor. This in turn, by Equation 4.20, 
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is dependent on the specific power of the reactor and varies 
inversely with the concentration of uranium in the core, so 
that low power density and high uranium concentration are 
desired . 

As is shown in Section 6.1, the operating power of 
the reactor is proportional to the area of the core, with 
core temperature and pressure fixed. The specific power varies 
then inversely as the core radius, as does the concentration 
of uranium in the core at small radii, shown in Figure 7.. 

It can be concluded that, at least, for many possible 
re-entry materials there would be no difficulty in handling 
the models after testing and that model contamination is not 
a critical parameter in the design of this system, but should 
a particular material to be tested become highly contaminated, 
the degree of contamination can be reduced by reducing either 
the test duration or the specific power density. It is assumed 
that this does not occur in this study. 
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CHAPTER V 



REACTOR COOLING 



5.1 Heat Loads 

The walls of the 'gaseous core reactor are subjected to 
high heat loads due to convective, conductive and radiation 



solid materials of the core materials. 

At the high temperatures of the core the heat load due 
to radiation transfer is much greater than the convective 
and conductive heat loads: the convective and conductive 
heat loads are considered negligible in this analysis. The 
core is assumed to radiate as a black body at 10,000°K and 
the walls absorptivity set as 0.6. 

The heat load due to nuclear reactions is assumed to 

be due to neutron and gamma fluxes from the core and is set as 

9 

eight percent of the power output divided by the core area . 
The heat load per unit wall area is 



heat transfer and heat producing nuclear reactions in the 




(5.1) 



where 



6 is the wall absorptivity, and 
<T* is the Stephan-Boltzmann constant 
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For the conditions of this study, the heat load to the 
wall, given by Equation 5.1 is 

= S. 43 * lO ^Vcm 1 - sec. 

5.2 Criticality Effects 

The loads given by Equation 5.1 are much too high for 
the thick beryllium moderator to conduct away from the core, 
and are predominantly due to radiation. Reference 7 suggests 
the use of a thermal barrier, cooled by transpiration 
and radiation seeding to protect the walls of a gaseous core 
nuclear rocket under similar high radiation heat loads. 

This method of cooling is not appropriate for use in the 
reactor used in this system. The large size of the reactor 
is due to the high absorption cross section of the working 
gas, nitrogen . This causes a low value of the infinite 
multiplication factor in the reactor, Kc© • given by 
Equation 3.9 . Examination of Equation 3.9 shows that 
increasing the macroscopic absorption cross section in the 
core will reduce the infinite multiplication factor, . 

As this is already very close to 1.0 in this reactor, addition 
of even a small thermal barrier causes the reactor to be 
subcritical for any reactor size. 

Even if a material, with sufficiently low neutron 
absorption characteristics and high strength could be found, 
the transpiration flow necessary to internally cool the 
thermal barrier is not available due to the large area of 
this core and the relatively small throat area. 
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The moderator in this study is protected from the 
high thermal radiation heat loads by a seeded layer of carbon 
particles in the core near the walls. The small amount of 
‘carbon necessary to protect the walls does not preceptibly 
lower the multiplication factor, but has a more pronounced 
effect on lowering the diffusion length, defined by Equation 
3.8 , so that the carbon acts as a moderator and slightly 

decreases the critical size, as shown by Equation 3.10 
This effect is due to the very small absorption cross section 
of carbon, which effects k» , and the larger scattering 
cross section, which reduces Ij . The effect of the carbon 
seed on criticality is ignored in this analysis. 

5.3 Reactor Flow Pattern 



In order to protect the moderator of the core, it is 
necessary to establish the following flow pattern. Consider 
any diagonal plane which includes the axis of symmetry of 
the nozzle. The reactor core flow and the nozzle are axi - 
symmetric. The reference direction,© - O , is aligned with 
the nozzle axis; the flow has a source at e-tr , the fuel 
and nitrogen inlet, and a sink at the nozzle. A boundary 
layer -is formed along the walls as the flow proceeds from 
inlet to the nozzle, much like a low speed divergent-con- 
vergent nozzle. 

The uranium fuel, in the form of a fine powder mixed 
with nitrogen, is sprayed into the core at the 0 ■=■ TT’ 
position. The nitrogen is sprayed into the core through a 
annulus around the fuel nozzle. Perfect mixing of the fuel 
and nitrogen is assumed. It is assumed that the carbon seed 
tends to remain toward the walls and is limited to a layer 
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one-tenth of the core radius thick and starting at the core 
wall. 



The gamma and neutron fluxes to the moderator is un- 
affected by the thermal radiation barrier. The total heat 
load to the moderator is the sum of the nuclear heat load 
and the thermal radiation which is not absorbed by the seeded 
layer. The moderator is cooled by a combination of radial 
injection of a fraction of nitrogen flow into the core, and 
water cooling. The nitrogen flow itself is not sufficient 
to cool the moderator and must also be used to establish the 
desired flow pattern. 

5.4 Radiation Barrier 



With the flow pattern in the core as established in 
paragraph 5.3, the flow is essentially flow in a nozzle 
with a diverging section proceeding the converging-diverging 
nozzle. Reference 12 gives the results of using a seeded 
layer to protect a gaseous core rocket nozzle. With perfect 
mixing of the fuel and working fluid, a seeded layer 1/10 
of the local nozzle radius thick will reduce the radiant 
energy flux by two orders of magnitude if the seeded layer 
absorption coefficient. It , is 1/3 centimeter.^" 

Reference 13 gives a correlation of the total absorption 
coefficient per seed particle concentration, K / ^ , and the 
particle radius. Using carbon particles as the seeding 
material, it is desired that the average concentration of 
carbon particles be minimized, both for nuclear and thermo- 
dynamic reasons. 
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The average carbon density, 0 , in the seeded layer is 



where 




(5.2) 



is the mass density of the carbon seed 
is the average seed radius 

With K set at Y& Cm / the average density is minimized 

at 3.7 1 10 7 ram /centimeter 3 *, this occurs with a concentration 
»© 3 

of 4.'5 •* JO particles/cm having an average diameter of 

0.1 microns. 



ft 

and R 3 



As the seeded layer does not fill the core, the average 
seed density over the core is 



fc - 



fiir (r?-9R I) 



(5.3) 



tr R? 



a i<IO T ra.ms /ce« tlmefcer* 



The contribution to the macroscopic absorption cross 
section due to the carbon atoms in the core is given by 



£ 



?c cA 



SL 



Ji c < 5 - 4 > 

where cJL is Avogadros number and A c is the molecular weight 
of carbon* 



For the concentration of carbon in the core due to the 
seeded thermal barrier. Equation 5.4yields an absorption 
cross section contribution of \.G> * \0 7 cenU moVec”' . As 

this is four orders of magnitude less than the absorption 
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cross section of the nitrogen and uranium mixture in the 
core, the presence of the .carbon seed on criticality is 
ignored . 

5.5 Moderator Cooling 

The seeded radiation barrier discussed in Section 5.4 
reduces the radiant heat load to the wall, given in Section 
5.1, by two orders of magnitude, but does not effect the nu- 
clear radiation induced heat load. The resulting heat load 

2 

to the moderator is 350 cal/cm sec.. 

The thickness of the moderator and the thermal conductivity 
of beryllium prevent this heat load from being conducted 
through the moderator and removed by back side cooling. 

A portion of the incoming nitrogen flow is radially 
passed through the moderator to absorb the heat load to the 
moderator. The coolant nitrogen flow is injected tangentially 
into the core in such a manner as to be parallel to the flow 
of the boundary layer at the point of injection. This is 
to prevent the injected flow from driving the seeded thermal 
barrier away from the wall. A small amount of carbon, as 
is necessary to replace carbon in the thermal barrier which 
vaporizes or diffuses out of the thermal barrier, is carried 
by the cooling nitrogen flow. 

The maximum temperature of the beryllium moderator is 
set at 1000°K by structural strength limits -With one half 
of the primary nitrogen flow used to establish the desired 
flow pattern, the remaining primary nitrogen flow is capable 
of absorbing one third of heat load in the moderator, increasing 



47 



. o o . 

in temperature from 293 K to 1000 K- The remaining heat 

load is absorbed by bleeding four percent of the water flow 

enroute to the water diffusor and allowing its temperature 

o 

to increase to 1000 K in the moderator. 
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CHAPTER VI 

NOZZLE HEAT TRANSFER 



6.1 Nozzle Size 



The hot gases in the core are radiating energy away 
from the core. As the moderator walls are not capable of 
transferring this heat load, they must be protected by a 
seeded thermal barrier, and the absorbed energy carried out 
of the core through the nozzle. Thus the minimum operating 
power is set by the radiant heat load from the core, and the 
size of the core wall surface. The nozzle throat area is 
found by equating the radiative power in the core to the 
energy flow through the nozzle, then. 



A* = P-2^_Ai 

* f-V-Vlo 



( 6 . 1 ) 



For the operating conditions set in Section 2.2, the 

. 2 
minimum throat size is 93 centimeters . 



6.2 Heat Loads 



Large heat loads are experienced by the walls of the 
nozzle as the flow is accelerated from stagnation conditions 
of the core to the high Mach number desired in the test 
section. The heat load is due to convective and radiative 
heat transfer from the hot gases and nuclear heating of the 
nozzle walls by alpha, beta and gamma interaction. To 



49 

prevent failure of the nozzle, this heat load must be 
transferred through, absorbed by or blocked from the nozzle 
walls. 

Turbulent flow is assumed in the nozzle, also it is 
assumed that the flow is frozen, isentropic and one dimensional 
in determining local flow properties. The heat generated in 
the nozzle by nuclear interaction is assumed negligible compared 
to the convective and radiation heat loads. 

6.3. Nozzle Configuration 

The nozzle consists of a series of longitudinal porous 
tubes, the area ratio of which is adjusted to provide an 
over pressure, of the cooling fluid in the nozzle tubes 
above the pressure of the flow in the nozzle, of 10 atmospheres. 
This provides the driving potential for transpiration cooling 
of the nozzle. Flow inside the coolant tubes is assumed to 
be turbulent. 

The tube cross section is a "U" shape, the flat sides 
of which may vary to provide the desired area ratio and also to 
facilitate joining the tubes together to form a smooth nozzle 
shape. The heat transfer surface is approximated by a circle 
whose radius is R. The fin effect of the radial sides of 
the cooling tubes is not accounted for, giving a conservative 
value for the heat transfer to the coolant. 

The coolant tubes support the thermal stress due to 
heat transfer across their radially inner-most face. They 
also support the hoop stress due to the coolant pressure. 

The tubes are closed by a shell which supports the pressure 
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of the hot gas in the nozzle. 

6.4 Maxirr.ura Heat Transfer 

The maximum heat load that can be transferred through 
the wall is set by the material properties, the configuration 
of the wall and the pressure the wall must support. 

For a thin walled tube, with' the wall thickness, d, very 
much less than the tube radius, , the temperature gradient 
across the tube is linear. The thermal stress across the 
tube, in the tangential direction, are maximum in tension at 
the cool inner fiber, 

( T w -Tw* )/e (6.2a) 

and maximum in compression at the hot outer fiber 

= _< ^L ( T ^-Tw 0 ^/2. (6.2b) 

»vg 

where 

C( is the coefficient of linear expansion, 

E is the modulus of elasticity, 

is the hot side wall temperature, 

T^is the cold side wall temperature, 

and both 0C and E are functions of temperature. 

■ The pressure differential across the cooling tube wall, 
■p, causes a uniform hoop stress in the wall, 

* = t ±3x. 

A (6.3) 

where 

d. is the wall thickness. 
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and the positive sign indicates tensile hoop stress, the 
negative compressive. 



The total stress is given by the sum of the thermal 
and hoop stress; scaler sum in this case as both are tangential. 



At the extreme fibers, 

<r= - qcE'CTL-Xio) v Rt 

2 a 



(6.4) 



The configuration chosen in this analysis utilizes a series 
of small cooling tubes forming the nozzle. This gives a 
tensile hoop stress and the maximum total stress is tensile at 
the cool inner fiber. The hoop stress relieves the compressive 
thermal stress at the hotter inner wall, where the material is 
weakest . 



The heat transferred through the wall is given by. 



where 



0 = y (Tv.-'Lo) 



k is the wall thermal conductivity 



(6.5) 



Equation 6.4 is solved for the temperature across the 
wall, and this is substituted into Equation 6.5 , 




- -PRt 
[ ~ ~5~ 

KE d. 

2k 



( 6 . 6 ) 



The maximum heat transfer through the wall, with respect to 
wall thickness, d, is obtained by setting, dQ /dcL equal zero 



4kiP Rt 

oCEcP 



2<rk _ 

a£d l 



O 



(6.7) 
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The condition for maximum heat transfer through the 
wall is. 



5k 

2 ' dL 






( 6 . 8 ) 



The maximum heat flow through the wall is 

Q-*** ~ 2^ ET f * 6,9) 

Here failure bccurs when the tensile stress at the inner 
fiber of the coolant tubes reaches the ultimate stress of the 
tube material. Radial stress is due to the pressure differen- 
tial across the tube thickness and is negligible compared to the 
tangential stress. The nozzle has incroporated an expansion- 
slip joint at the up stream, low Mach end; this prevents 
longitudinal thermal stress so that the stress in the coolant 
tubes is considered uniaxial. 



6.5 Material Selection 



A material was sought with a maximum value of *"*g= 

in order to maximize heat transfer with respect to material 

selection. This parameter is a function of temperature in a 

given material and the temperature must be high enough so 

that super-cooling is not necessary at the back side to absorb 

the heat load. Wall temperature also affects the heat load 

10 

to the wall but examination of the Bartz's Equation , 

utilized to determine the convective heat load, shows only 

a weak dependence of convective heat transfer on wall temperature. 

24 o 

This investigation chooses Rene 41 at T w of 1025 K 
as having the desired structural properties. 



53 



At this temperature, 

2 

Ci. is 146,000 pounds/inch 

_ . 7 ,2 

b is 2.4 x 10 pounds/inch 

K. is 0.05 cal/cm-sec-°K 
-7 , o 

1.57 x 10 cm/cm K 

Minimum wall thickness is set by manufacturing limitations 
and is assumed to be 0.01 cm for this analysis. 

This choice of material and wall thickness gives a 
maximum heat transfer per unit wall area of 

* 1-94 * flej./cm 1 sec 

provided that the cooling tube radius, set by Equation 
6.8 with a ten atmosphere pressure differential across the 
wall, is 5 cm or less. The temperature drop through the 
wall is 

T. - K,. - -589”\c 

and T w „ = C3C*k. 

6.6 Test Section Mach Number 



The nozzle shape considered is a convergent-divergent 
nozzle formed by the previously defined cooling tubes. 

The nozzle is axially symmetric with a 45 degree convergent 
section, 15 degree divergent section and a radius of curvature 
at the throat equal to the throat inner radius. 

The heat load to the nozzle is weakly dependent upon 
throat radius, affecting the convecting heat transfer co- 
efficient only in this analysis. The pump work required of 
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the system as a whole, however, varies directly with the 
mass flow, which varies as the area of the throat. 

The throat radius is chosen in Section 6.1 as 5.45 
centimeters . 



The test section Mach number is arbitrary with a friction 
less isentropic nozzle. In this nozzle, transpiration cooling 
is required until the flow has expanded to about Mach 6, 
which requires a lengthy nozzle, of the order at 15 meters. 
Expanding the flow to Mach 6 increases the transpiration 
coolant requirement and nozzle losses. For this analysis, a 
test section Mach number of 3.0 is chosen to enable simulation 
of the desired flow fields, but reduce transpiration pumping 
requirements and nozzle losses. 

6.7 Transpiration Cooling Requirements 



The heat load per unit area to the wall, neglecting 
gamma heating, is given by, 

Q - (6.io) 

where Qo* is the convective heat load in the absence of 
transpiration cooling, 

Qr* is the radiative heat load, assumed unaffected 
by transpiration cooling. 

The convective heat load was obtained by the method of Bartz' 
for turbulent compressible flow. 



where 



Q©c - V\ (, Taw -TL ) 




(6.11a) 



9 

<r 



(6.11b) 
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£T= 



1 




Uj 

5 



(6.11c) 



V\ is the heat transfer coefficient, 

4J is 0.6, 

D is the throat diameter, inches, 

J£ is the viscosity, 

Cp is the specific heat at constant pressure 
nr is the temperature, °R , 



sub 



wall 



w 



sub 



bulk 



o 



sub 



adiabatic wall 



aw 



is the chamber pressure, 

t*. is the radius of curvature of throat, inches. 

It is 1 for the nozzle, 

is the nozzle area ratio, 
y is the ratio of specific heats, 

M is the Mach number. 

14 

The thermodynamic properties of nitrogen were used in 
this analysis, and the presence of the uranium was ignored. 

The radiation heat load is obtained by considering the 
gas to radiate at stream ambient temperature as a black body 
and the wall to have an absorptivity of 0.6, 



The radiation from the wall is negligible compared to the 
load to the wall and this reduces to 



Q* * c* u - e* <r O 



( 6 . 12 ) 



Q* * 6“ T A 



(6.13) 
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where 

is the wall absorptivity, 

€* is the gas absorptivity, 

CT* is the Stefan-Boltzmann constant, 

T is the Temperature, ambient free stream 
is the wall emissivity. 



The heat load to the nozzle wall in the absence of 
transpiration cooling is given in Figure 8 . The maximum 

heat load that the wall is capable of withstanding is lower 
than the load experienced. Transpiration cooling by the 
injection of low temperature nitrogen is used to reduce the 
convective heat load to the wall. 



The effects of transpiration cooling in reducing the 
convective component of heat transfer is given in Referencell. 
The results are given in terms of the Stanton number reduction 
obtained by injection of a homogeneous coolant into a turbulent 
nozzle flow verses a blowing parameter , F/st®. Here F is de- 
fined as the ratio of the transpirant mass velocity to the 
local free stream mass velocity. 



F = 



(6.14) 



where G is defined as the transpiration mass velocity. The 
Stanton number is defined as. 



h (6.15) 

?V c e 

Then the Stanton number reduction is, to the first order, 
the reduction in convective heat load. 




St / Sr 0 — Q.c / Qco 



.(6.16) 
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The maximum heat load transferred through the wall, 
with transpiration cooling, is the sum of the reduced convective 
heat load and the radiation heat load minus the heat absorbed 
by the coolant passing through the wall. 

CU. = Q, » Q*- Tw'-O (6-17) 

where Cp is evaluated at local pressure and the average 
temperature in the wall. 



The convective heat load without transpiration cooling 
is given by, 

Qc. = U. (T^ -X.) (6.18) 



and 



Me. , -EL 

ho 



(6.19) 



With Equations 6.18 and 6.19 the maximum heat flow through 

the wall can be rearranged to 

Qwm _ Q » _ Oc ^ ( i^«* o q \ 

Qoc Q 0fc ' Qo. s x . c p (T^-T^ V ' 

The left hand side of this equation is known at each 
station and the right hand side suggests the possibility of 
a graphical solution using the difference between the Stanton 
number reduction data of Reference 11 and a family of straight 
lines of slope C p (T** -7 W .) / C ? ( Taw -T w ) . Figure 9 

reproduces the Stanton number reduction data of Reference 11 
and illustrates the graphic method of solution. The blowing 
parameter, F, is obtained by this method and from it, the 
transpiration mass velocity necessary to reduce the heat 
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load through the wall to the maximum allowable. This flow 
is graphically integrated to obtain the total required 
transpiration flow. 

The transpiration cooling flow per unit axial nozzle 
length, 2ttRGi / is given in Figure 10 , along with the total 

transpiration mass flow. 

Under the conditions of this study, the transpiration 
coolant flow necessary to prevent thermal damage to the 
nozzle and test section is l.25x \0* ya.m/scco«<i , and 
the ratio of transpiration flow to the primary flow in the 
nozzle is, 

= 0.37 

This estimate on transpiration cooling requirement is 
conservatively high as the flow is assumed to radiate as a 
black body and the effects of thermal radiation blockage by the 
seeded layer in the core is neglected. 

6.8 Secondary Nitrogen Flow ' 

The heat flow through the wall must be absorbed by 
the flow in the coolant tubes. The flow properties in the 
tubes are set by matching the heat transfer at the throat 
as this is the point of maximum heating load. 

The flow in the tubes can not be set point by point 
to enable maximum heat transfer to be possible along the 
entire length of the nozzle and heat transfer rates below 
are possible along the nozzle. This problem is 
recognized but not investigated further in this analysis 
as it is felt that the transpiration cooling flow can be 
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tailored to fit a wall with variations in heat transfer 
capabilities . 

The heat flow through the walls is absorbed by the 
coolant flow, 

GU. * Vu (Tw-Tc) (6.21) 

Here radiation and nuclear heating of the coolant are ignored. 
The coolant is flowing in a nozzle shaped duct but here 
fluid is bled out of the tube. In the absence of data on 
sucking of a turbulent boundary layer, the data of Reference 11 
is extrapolated to include the sucking case. This is also 
illustrated in Figure 10, and is analogous to the treatment 
of a sucked laminar boundary layer given in Reference 25. 

The method of Bartz is used as before to determine the 
heat transfer coefficient. The coolant tube nozzle throat 
radius is set at .5 cm. 

The chamber temperature and pressure of the secondary 
nitrogen flow is set by matching the heat transfer through 
the wall at the throat to the heat absorbed by the coolant 
of that point, taking into account the effect of the transla- 
tion cooling flow sucking the boundary layer in the coolant 
tubes. 

* 

For the conditions of this study, the secondary nitrogen 

chamber pressure is 1010 atmospheres and the temperature is 
o 

530 K. With these chamber conditions the mass flow of 

5 

secondary nitrogen is 2.17 x 10 grams per second and the 
ratio of secondary to primary nitrogen flow is 0.67 to 1. 
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The secondary flow of nitrogen that is used to provide 
transpiration and backside cooling of the nozzle is approx- 
imately the same as the primary flow through the nozzle. This 
doubles the amount of pump work required, but as the nozzle 
would fail without the transpiration protection, this additional 
pump work is an acceptable penalty that must be incorporated 
into the design of the system. 
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CHAPTER VII 
AUXILIARY EQUIPMENT 

7.1 General 

Operating this facility as a closed cycle removes 
some of the technological problems involved in the design 
of a gaseous core rocket. The problem of containing the 
gaseous fuel in the reactor core is avoided here by operating 
as a closed cycle. This introduces the problem of additional 
equipment in the cycle. Among the additional equipment needed 
are high powered turbines, high capacity pumps, large heat 
exchangers and fuel separation and processing equipment. 

7.2 Water Diffusor 



The flow at the exit of the test section is mixed 

with the secondary nitrogen flow which was utilized to provide 

transpiration and backside cooling of the nozzle. At this 

point the flow stagnation temperature, determined from the 

o 

enthalpy of the mixture, is 8000 K and the stagnation pressure, 
assuming isentropic expansion through the nozzle, is 1000 
atmospheres .Conventional heat exchangers can not operate 
at these temperatures and pressures, neither can conventional 
turbo-machinery. To lower the stream's temperature and pressure, 
a water spray diffuser is employed . The water is sprayed 
directly into the hot nitrogen stream and is heated to super - 
heated steam by the nitrogen, which loses energy in the process. 
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The cross sectional area of the diffuser can be designed 
to obtain the desired outlet conditions. 

Conservation of energy determines the required water 
flow. Conservation of energy requires that 






where yy\ is the mass flow rate 
H is the total enthalpy 
1 is. the mole fraction of water 

1’^ is the mole fraction of nitrogen 
and subscript w refers to water 
i inlet conditions 
0 stagnation conditions 
2, outlet conditions. 



(7.1) 



Equation 7.1 was utilized to obtain the required mass 
flow of water necessary to reduce the temperature and pressure 
of the nitrogen stream to the desired values. 



The size of the water diffuser is set by the perfect 

gas law, written in the form 

* 

P _ A. yf\ i [ t a y, o'K i (7.2) 

P 2 A z T, V x cH. x 

where all values are static values and 

Ji is the mixture average molecular weight 

^ is the mixture average specific heat ratio 

With the properties of nitrogen given in Reference 
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14 and those of water from Reference 26 , Equations 7.1 
and 7.2 can be used to find the necessary water flow and 
diffusor exit area, once the exit Mach number, temperature 
and pressure are specified. 

The diffusor outlet total temperature is set at 1000°K, 
as an upper operating temperature for turbo-machinery. The 
outlet total pressure is set a 100 atmospheres and the 
Mach number at 0.01 

With these conditions set, the required mass flow of 
5 

water is 9 x 10 grams per second, and the output area 

6 2 

of the diffusor is 8.25 x 10 cm. 

7.3 Fuel Separator 

The cooling process in the water diffuser cools uranium 
fuel in the flow well below its melting point. The super- 
cooled uranium rapidly condenses and forms small droplets 
and particles which are suspended in the flow.® The thermal 
energy of the flow is used to provide a centrifugal force 

field to separate the particles from the gas. The drag on 

25 

these particles is given for Stokes flow as 

£> = **>L V (7.3) 

m 

where is the viscosity of the steam, nitrogen mixture. 

The equation of radial motion of a uranium particle in 
a centrifugal force field is, with the multiplier of velocity 
in Equation 7.3defined as Co, 

r + C© f - os* r = O (7.4) 



This has a solution of the form 
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') 



(7.5) 



For small particles, the drag coefficient C© is 
large and the second exponential term of Equation 7.5 
decays rapidly regardless of the value of <x) and the 
first term determines the radial position of the 
particle. For small values of co , the particle is strongly 
restrained by the drag forces of the fluid and a centrifugal 
separator is not very effective. For a centifugal separator 
to efficiently separate the small particles, the rotational 
velocity must be at least of the order of the drag cofficient. 
Under the conditions of the flow just after the water diffusor, 
the angular velocity must be of the order of 100 radians 
per second for micron sized particles. This imposes large 
velocity requirements if the flow is swirled in one section 
at the outlet of the water diffusor, so the flow is divided 
into several centrifugal separators. 

7 .4 Turbines and Pumps 

The ideal power available in the flow as it comes out 
of the uranium separators is given by 



It is convenient however to expand the flow only until 
the water vapor in the flow reaches its saturation point. 
At this point, the water can be condensed and removed from 
the flow stream and the nitrogen further cooled to reduce 
the pump work required to increase its pressure to core 
entry conditions. Under these conditions Equation 7.6 



JP - m* H w 



(7.6) 
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becomes 

IP - mM |^14 n (lOOo’K/t-SJiOOa.Vm.) 

+ W*, ('OOOK, i\OOaV«,) 

For isentropic expansion of the water vapor, the ideal 
power available is 2,520 megawatts. 

After condensation, the water is pumped as a liquid 
to 120 atmospheres and injected into the water diffusor. 
This requires an ideal pump work, 

O 

■W Pw - (7.8) 

Ideally, 10.9 megawatts are required for the water pump 
work. 



- U-(WT f U'l] 

- P.„) J (7.7) 



The nitrogen is cooled in a heat exchanger to 293°K 
and then pumped back to 1000 atmospheres. This must be 
done in at least two pumps with interpump cooling to maintain 
the temperature at or below a structurally safe limit. The 
ideal nitrogen pump work is given for each pump by 

= [vU\,PO-H*(T t) R.)] (7.9) 

where the temperatures and pressures are related by the 
isentropic compression relation 



v-t 




(7.10) 
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The total pump work for all nitrogen pumps is 868 
megawatts . 

The total heat rejected in the cycle is the sum 
of the heat of condensation of the water vapor, the 
cooling of the nitrogen prior to being pumped, and inter- 
pump cooling after pumping so that the nitrogen can be 

used to cool the reactor moderator. The total heat rejected 

3 

in the cycle is 7.1 x 10 megawatts. 



If water is used as a heat sink for this rejected heat, 
the mass flow required is 






Q 






dp rTT 



(7.11) 



o 

For a temperature rise of this coolant water of 100 K, 

g 

the mass flow required is 1.1 x 10 grams per second. 



7.5 Minimum Power Level 



The minimum operating power level of this test facility 
is set by the criticality requirements of the reactor. Chapter 
III, and heat transfer considerations in the reactor and nozzle. 
Chapter V and VI. This high operating power requires high 
mass flow rates and large capacity pumps and turbines, as 
indicated in Section 7.4. 

I • 

There is no lack of power available to drive this machinery 
and small units can operate in parallel to provide the desired 
capacity. As the system is proposed, there is a requirement 
that the combined efficiency of the turbines and pumps be- 
about 35 percent, but if all of the energy in the cycle is 
used, efficiencies as low as 12 percent are acceptable. 
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Equations 5.1 and 6.1 can be combined to express the 
operating power minimized "with respect to the area of the 
throat , 

£> = (€G-T 0 4 + .08 P/AJ Ac (7.12) 

The area of the core is minimized for any pressure 
and temperature by setting the ratio of nitrogen mass 
flow to uranium mass flow at 15.4, as was shown in Section 3.2. 

Equation 3.11 shows that, with the mass flow ratio 
fixed, the critical radius, R, varies inversely with the 
number density of uranium, N„ . For constant operating 
pressure, the number density of uranium varies inversely 
with the operating temperature, 1^, . The critical radius 

is then proportional to the operating temperature and the 
core area is proportional to the square of the operating 
temperature . 

With fixed mass flow ratio and operating pressure, the 
minimum operating power is, to the first order, 

P„.„ = C T 0 ‘ (7-13) 

where C is a proportionality constant. 

Equation 7.13 is valid only in the temperature range 
where the heat loads to the moderator are too large to be 
removed by conventional cooling and the energy absorbed in 
protecting the moderator walls must pass through the nozzle. 



A plot of the estimated minimum operating power as a 
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function of temperature is given in Figure 11. At operating 
temperatures below 7000°K, this facility is in the power 
category of large arc tunnels and MHD devices, and the 
relative merits of each type of facility would have to be 
considered before deciding which system is most desirable. 

At the higher temperatures at which this facility is con- 
ceptionally most useful, the minimum operating power rises 
rapidly, as well the size of the physical plant. It is also 
noted that reducing the operating temperature reduces the 
re-entry simulation capability of the test facility, as 
illustrated in Figure 2. 

The economics of developing such a facility must be 
examined to determine if the need for such a testing capability 
justifies the expense of so large an undertaking. 
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CHAPTER VIII 
CONCLUSIONS 

Within the scope of the material examined in this 
report it appears to be within the capability of projected 
technology to design and build a high enthalpy test facility 
powered by a gaseous core reactor. 

This test facility can not exactly duplicate all of 
the conditions imposed upon a vehicle during an atmospheric 
re-entry, but does provide a means of testing re-entry 
materials and shapes under the high temperature and shear 
conditions which are encountered during a re-entry. Existing 
facilities cannot attain the high temperature, high pressure, 
long duration flows necessary to test the shear mode of 
failure. The gaseous core reactor provides the power required 
for these flows. 

Nuclear contamination of several proposed re-entry 
materials is slight enough so that the model may be examined 
immediately after a test run. This may not be true in 
general and each proposed test material must be examined 
by the methods of Chapter IV to determine the degree of 
contamination . 

The high neutron absorption cross section of nitrogen 
as a working fluid causes the reactor to be large and the 
mass flow ratio of uranium to nitrogen to be relatively 
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high. The reactor size and operating power is minimized 
by setting the mass ratio of nitrogen to uranium in the 
core at 15.4 to 1. This adversely affects the quality 
of the mixture for simulating atmospheric conditions, 
but is slight enough to be considered acceptable. 

The minimum operating power of the test facility is 
set by the radiant power of the core which cannot be 
absorbed by the moderator and must pass through the nozzle. 

In the neighborhood of 10,000 °K, for a constant mass ratio, 
the minimum operating power varies as the temperature in 
the core to the sixth power and as the inverse of the core 
pressure . 

Operation of the facility at the minimum power level 
requires high mass flow rates and associated high capacity 
pumps and turbines. These pumps must be capable of high 
pressure and high pressure ratio operation. Efficiency 
of the combined turbine and pump can be as low as twelve 
percent for self sustaining operation of the system. 

The turbo-machinery required set a practical upper 
limit on the operating power of the system. Structural 
considerations impose an upper pressure limit on the system. 

This analysis has examined only the minimum powered 
point for the operation at the maximum pressure and temperature 
and has found that such a facility appears technically feasible 
Operation of a facility at lower power levels is possible 
with a reduction in chamber temperature and associated loss 
of re-entry simulation capability. 



71 



It is recommended that the economics ' involved in the 
construction, maintenance and operation of this system be 
examined to determine if such a large facility would be 
economically practical. 
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TABLE I 

FISSION PRODUCTS 



Radionuclide 


Half Life 


Yield Fraction 


Decay Constant 


A 




t»/a (s*c ) 


*1- 


Aa ( see"' ) 


Br 


87 


55.6 


10- 3 


1.2 x 10" 2 


Br 


89 


4.51 


5 x 1CT 3 


1.54 x 10 _1 


Ag 


110 


24.2 


10- 4 


2.9 x 10" 2 


In 


114 


72 


io " 4 


9.6 x 10 -3 


In 


116 


13 


10 - 4 


J -2 

5.4 x 10 


Sb 


135 


1.52 


4 x 10~ 2 


4.56 x 10 _1 


I 


137 


22 


6 x 10' 2 


3.1 x 10" 2 
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STAGNATION ENTHALPY 
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FIGURE I REENTRY STAGNATION 
PRESSURE AND 
ENTHALPY. 



POWER PER UNIT THROAT AREA, 7a. , ME6AWATTS/CM 
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INITIAL ENTRY VELOCITY, V e , ( I0 3 M/SEC) 

FIGURE 2 REENTRY SIMULATION 
POWER REQUIREMENTS. 
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FIGURE 3 TEST FACILITY SCHEMATIC. 
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4-5 
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n 2 pump 


7-8 


COOLER 
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n 2 PUMP 


9-10 


COOLER 


10- 1 
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FIGURE 4 NITROGEN CYCLE. 



TEMPERATURE (10^° K) 
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FIGURE 5 WATER CYCLE. 



-78- 




( £ W0/Sl/\101V) N0I1V81N3DN00 IAIDIN ViJD IVOIlldO 



4 



FIGURE 6 CRITICALITY OF BARE URANIUM-NITROGEN GASEOUS REACTOR. 
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FIGURE 7 CRITICAL RADIUS OF BERYLLIUM REFLECTED URANIUM- NITROGEN 
GASEOUS REACTOR. 
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FIGURE 8 NOZZLE HEAT LOADS. 
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FIGURE 9 TRANSPIRATION COOLING 
REQUIREMENT. 







FIGURE 10 NOZZLE TRANSPIRATION 
FLOW. 



ESTIMATED MINIMUM OPERATING POWER (MEGAWATTS) 
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FIGURE 1 1 MINIMUM OPERATING POWER FOR 
GASEOUS CORE REACTOR. 
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